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Abstract Results from a GPS reflectometry experiment with a 12 channel 
ground-based GPS receiver above two lakes in the Bavarian Alps are presented. 
The receiver measures in open- loop mode the coarse/aquisition code correlation 
function of the direct and the reflected signal of one GPS satellite simultaneously. 
The interference between the coherently reflected signal and a model signal, 
which is phase-locked to the direct signal, causes variations in the amplitude 
of the in-phase and quad-phase components of the correlation sums. From 
these amplitude variations the relative altimetric height is determined within a 
precision of 2 cm. 

1 Introduction 

Satellite-based active altimeters on ENVISAT and JASON deliver valuable 
ocean height data sets for global climate modelling. In order to improve the 
climate models, altimetric data of higher resolution in space and time is re- 
quired. This gap can potentially be filled with GPS-based altimetric measure- 
ments. Additionally, ground-based GPS receivers can monitor ocean heights in 
coastal areas wher e sate l lite altimetry data get coarse and decrease in quality 
l^u and Cazcnavtl lioOlt IShum et all 1199/1) . 

Since GPS altime try has been proposed as a novel remote sensing capability 
ijMartfn-Neiral Il993l) . many studies have been carried out at different obser- 
vation heights and platforms. While Earth-reflected GPS signals have been 
observed from spaceborne instruments ijLowe et al.L l2002at iBeverle et al.l l2002ll 
and the CHAMP and SAC-C satellites already are equipped with Earth/nadir 
looking GPS an tennas, work is in progr ess in order to establish satellite-based 
GPS altimet rv llHaii and Z uffada, 200^ 1). Airborne campaigns have been con- 
ducted (e.g. iGarrison et al.. ((1998.) . .Garrison and Katzber Jll2f)0(Tl). iRius et alJ 
l|2002l) ) and recently reached a 5-cm height precision l|Lowe et alll2002bl) . Ground- 
based GPS altimetry measurements have been p erformed at laboratory sca le of 
some meters height with 1- cm height precisi o n dM artin-Neira et al .l l2002l) up 
to low-altitudes height (e.g. lAndersonI ll2000l). iMartm-Neira et all l)200l| Vr and 
reached a 2-cm height precision (jTreuhaft et al.U200l|) . 
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In this study a 12 channel GPS receiver is used llKellev et all 120021) . that 
has been extended with a coarse/acquisition (C/A) code correlation function 
tracking mode. In this coherent delay mapping (CDM) mode the direct GPS 
signal is tracked while concurrently the reflected signal is registered in open-loop 
mode. Using the LI carrier phase the relative altimetric height is determined 
from the components of the reflected signal. 

2 Experimental Setup and Data Acquisition 

The experiment was conducted on 8 - 10 July 2003, 50 km south of Munich, 
Germany, in the Bavarian alpine upland at the mountain top of Fahrenberg 
(47.61°N, 11.32°E) at a height of about 1625 m asl. Mount Fahrenberg belongs 
to the Karvendel mountains and from the mountain top unobstructed view is 
available to lake Kochelsee (surface area about 6 km^) to the north and lake 
Walchensee (surface area about 16 km^) to the south. Following a schedule of 
predicted GPS reflection events, the receiver antenna was turned towards the 
lake surface of Kochelsee (about 599 m asl) or Walchensee (about 801 m asl). 
The antenna was tilted about 45° towards the horizon. During a GPS reflection 
event the direct and the refl ected signals interfere at the antenna center (e.g. 
IParkinson and Spilkeil l)l99(j) ). The interference causes amplitude fluctuations 
that are quantitatively analyzed to determine the height variation of the specular 
reflection point. 

The receiver is based on the OpenSource GPS design llKellev et all l2002t) 
and was modified to allow for open-loop tracking of reflected signals. The receiv- 
ing antenna is an active RHCP patch antenna (AT575-70C from AeroAntenna 
Technology Inc.) with 4 dBic gain, 54 mm in diameter and a hemispheric 
field-of-view. Operating in CDM mode all 12 correlator channels are tuned to 
the same GPS satellite by setting the appropriate pseudo-random noise (PRN) 
value. The correlation between the r eceived and mode l (replica) signal is re- 
ahzed in hardware (Zarlink GP2021. EaRLINKI l|200lh 'l. While one channel 
(the master channel) continues to track the direct signal, the carrier and code 
tracking loops of the 11 remaining channels (slave channels) are synchronized to 
the master channel. Each channel consists of the prompt and the early tracking 
arm at zero and at -0.5 chip code delay, respectively. Thus, 2 x 11 = 22 delays 
are available to map the C/A code correlation function of the reflection signa- 
ture. In CDM mode the slave carrier and code phase-locked loops (PLLs) are 
opened and their loop feed-back is obtained from the master PLL. All carrier 
loops operate with zero carrier phase offset with respect to the master chan- 
nel; in the code loop feed-back, however, delays covering an interval of 2 chips 
(about 2 ^s) with a step size of 0.1 chips (about 100 ns) are inserted. In-phase 
and quad-phase correlation-sums of each channel are summed over 20 ms taking 
into account the navigation bit boundaries and stored together with code and 
carrier phases to hard disk at a rate of 50 Hz. Figure 13 illustrates the CDM 
mode: while the direct GPS signal is tracked with the prompt and early arm of 
the master channel at and -0.5 chips code offset, the prompt and early arms of 
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the remaining 11 slave channels are set to chip code offsets between 0.4 and 2.7 
to map the reflected signal (corresponding to an optical path difference of 120 
to 810 m). In Figure H the root -sum-squared in-phase and quad-phase values 
of the reflected signal are plotted as a function of code delay. The maximum 
power of the reflected signal is about 20 log 0.2 = —14 dB below the direct sig- 
nal's power. The peak of the correlation function is separated by a delay of 1.5 
chips from the direct signal's correlation peak. 

Data analysis is performed in the following way: first, the code delay cor- 
responding to the maximum of the reflected waveform is determined. Second, 
all in-phase and quad-phase correlation sum values Ir and Qr are extracted 
from the raw data which lie within a certain delay interval (grey box in Figure 

around the maximum code delay. The navigation message is demodulated 
according to 

ir ^ sign(/<j) Ir 

Qr = sign(/d)(3r, (1) 

where Id denotes the in-phase value of the master channel. Figure|3lA shows the 
oscillations of Ir and Qr caused by the interference between the reflected and 
the replica GPS signal. The phasor Ir + iQr rotates by about -I- 0.5 Hz due to 
the decreasing path length difference between the direct and the reflected signal, 
since during this measurement the GPS satellite moved towards the horizont. 
Note the phase offset of 90° between Ir and Qr- The phase (j) (Fig. 13 B) is 
calculated from the four quadrant arctangent 

= atan2(Qr,/r) (2) 

and is unwrapped by adding ± 27r when the difference between consecutive val- 
ues exceeds tt, resulting in the accumulated phase (pa- The optical path length 
difference 6 between direct and reflected signal is calculated from the accumu- 
lated phase (pa and the LI carrier wavelength A^i = 0.1903 m at the observation 
time t 

m-^>^Li- (3) 

Starting with a height estimate II{to), the temporal evolution of the altimetric 
height variation h{t) — /i(io)i normal t o the tangent plane at the reflection point 
P, is calculated from ( Trcuhaft ct alj l200ll) 

^ S{t)-S{to) + 2h{to) sina(to) 
^ ' 2 sina(i) ^ ' 

h{to) = {H{to) + rE) cos — - rE (5) 

rE 



with the arclength s defined in Figure^ an Earth radius — 6371 km and 

TT 

a = e + - - 7 (6) 

e = Cep/i + Aetropo, (7) 
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assuming an infinite distance to the GPS tran smitter. €eph is calculated from 
the broadcast ephemeris data l)GPS SPSlll995|) . the correction Aetropo accounts 
for refraction caused by atmospheric refractivity. The tropospheric correction 
is derived from a geometric raytracing calculation using a refractivity profile 
obtained from meteorological analyse provided by the European Centre for 
Medium-Range Weather Forecasts. The position of t he specular re f lection point 
P as function of 7 (Figure^ is calculated following ' Martin-Neiral l)l993(l . 

Thus, the altimetric height change of the GPS receiver above the reflecting 
surface is determined from the carrier phase difference between the direct and 
reflected GPS signal (Figure 13 C) . 

3 Data Analysis and Discussion 

During all 3 days several reflection events were observed from both lake sur- 
faces with different GPS satellites at elevation angles between about 10° - 15°, 
indicated by a clearly visible waveform (see Figure O. Several outliers can be 
observed in the data records. Most likely they are caused by overheating of the 
hardware correlator chip [S. Esterhuizen, University of Colorado, personal com- 
munication, 2003]. In this study outliers are removed in the following way: a 
value is calculated by linear extrapolation from the last 3 values of Ir{t). If the 
difference between extrapolated and actual value exceeds a threshold (20000- 
22000) , the extrapolated value is taken. The same is applied to the Qr (t) data. 
Additionally cycle slips (sporadic height jumps of about X^i m in adjacent data 
points) can be observed in the optical path length difference 6{t). The distor- 
tion of the data by cycle slips could be minimized by applying the same method 
as above to 6{t). Continuous data segments without cycle slips are chosen for 
height change determination. The mean receiver height above the lake surface 
is not expected to change during the short analyzed time periods. From to- 
pographic maps (scale 1:25000, Bayerisches Landesvermessungsamt, 1987) the 
heights H{to) are estimated to be 1026 m ±5 m (Kochelsee) and 824 m ±5 m 
(Walchensee) , respectively. By minimization of the linear trend of h{t) — h{to) 
we obtain a H{to) of 1022.5 m (Kochelsee) and 827.5 m (Walchensee). 

Figure 21 A and B plot the relative height change between the receiver and 
the reflection point at the surface of lake Kochelsee. Both observations used 
the same PRN, but were taken on different days. The height varies within an 
interval of about ±5 cm with a standard deviation of about 3.1 and 2.6 cm. 
Figure^lC and D show the height changes between the receiver and the reflection 
point at the surface of lake Walchensee. Again both observations were taken 
on different days and used different PRNs. Compared to the Kochelsee data, 
the height varies within a height interval of about ±2.5 cm with a standard 
deviation of about 1.4 and 1.7 cm. 

The different height variations at both lakes can be explained by different 
local wind and wave height conditions. As lake Walchensee is completely sur- 
rounded by mountains, waves are mainly driven by local, thermal induced winds 
which mainly occur at noon. Lake Kochelsee is open to the north, so longer 
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lasting wind can build up waves on the lake surface. 



4 Conclusions and Outlook 

Open-loop tracking of the reflected signals allows the determination of the rel- 
ative altimetric height with 2-cm precision. Different height changes can be ob- 
served at Kochelsee and Walchensee which reflect the different wind and wave 
height conditions at the two lake s. The relationship between t h e obse r ved height 
changes and win d speed (e.g. ICaparrini and Martin-Nei ra! Il9 98[). | Lin et all 
il999), Komiath v et all (1200(11) . IZuffada et all ll20n.'^ . ICardellach et all ll2nnlr i 
will be subject of further studies. The present receiver implementation is lim- 
ited to the observation of one GPS satellite at a time. To fully use the potential 
of GPS reflections the receiver will be modified to keep track of several GPS 
reflections simultaneously. 

Our results suggest that open-loop tracking is possible with low-gain and 
wide field-of-view antennas, showing the potential of this method also for space- 
based measurements of GPS reflections. 
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List of Symbols 



R receiver position 

P specular reflection point position 

e elevation angle of the GPS satellite above local horizon plane at R 
<5 observed path diflference between direct and reflected signal path 

rg Earth radius 
H receiver height 

h height variations normal to tangential plane at P 
a angle of reflection above tangential plane at P 

7 angle between normal of tangential plane and local horizon plane at P 
s arc length from subreceiver point to specular reflection point P 
in-phase correlation sum of the direct data 
Ir in-phase correlation sum of the reflected data 
Qr quad-phase correlation sum of the reflected data 
Ir Ir demodulated from navigation message 
Qr Qr demodulated from navigation message 

(p phase 
<pa accumulated phase 
Xli LI carrier wavelength 

^eph elevation angle calculated from broadcast ephemeris data 
Aetr-opo tropospheric correction to elevation angle 
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Figure 1: Geometry used to express the observed path difference S in terms of 
the known receiver position R with height H and the GPS sateUite elevation 
angle e and the calculated position of the specular reflection point P. h denotes 
the height variations normal to the tangential plane at P. Note that e has to be 
corrected by Aetropo due to the bending effect caused by the Earth's troposphere. 
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Figure 2: Delay mapped waveform of a reflection event (PRN f6) at 1334: f 7 
UTC 8 July 2003, antenna oriented towards Kochelsee. The delay is given in 
relation to the maximum peak of the direct signal. Blue circles and red triangles 
indicate 2 measurements (0.5-second duration) starting 120 (blue) and 267 (red) 
seconds after the start of the measurement. In the second case (red) the 2-chip- 
wide interval of covered chip code offsets is centered at the maximum of the 
reflected signal. The points reveal the measured waveform of the direct and 
reflected correlation signal. The thin black triangle marks the theoretical C/A 
code correlation function of the direct signal. The grey box marks the maximum 
of the reflected signal. 
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Figure 3: Panel A shows the demodulated reflected in- and quad-phase data Ir 
(blue circles) and Qr (red triangles) (PRN 16, elevation from 11.04° to 10.99°), 
antenna oriented towards Kochelsee at 1334:17 UTC 8 July 2003, as a function 
of time since measurement start. With Eq. [21 the phase (Panel B) and 
from Eq. ^andElthe relative height h{t) — h{to) is calculated (Panel C), with 
H{to) = 1022.5 m. 
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Figure 4: The left panels show relative height measurements at Kochelsee (PRN 
16), starting at 1334:17 UTC 8 July 2003 (Panel A) and starting at 1327:17 UTC 
10 July 2003 (Panel B) as a function of time since the start of the observation. 
PRN 16 changed elevation from 11.0° to 10.4° (Panel A) and from 11.4° to 
11.3° (Panel B). On the right panels height measurements at Walchensee are 
shown, PRN 20 (elevation from 14.7° to 14.1°), starting at 1257:15 UTC 8 July 
2003 (Panel C) and PRN 11 (elevation from 14.3° to 13.6°), starting at 1110:21 
UT 9 July 2003 (Panel D). H{to) = 1022.5 m (Kochelsee) and H{to) = 827.5 m 
(Walchensee). 
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